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RESUME 


Nous  presentons,  dans  ce  rapport,  une  thSorie  de  diffusion 
multiple  de  la  lumiere  pour  des  petits  angles,  derivee  de  1 ' interaction 
des  aerosols  sur  le  rayon  laser.  Le  calcul  de  la  luminance  Snergetique 
diffusee  peut  §tre  effectuS  a n'importe  lequel  degrS  de  multiplicity  de 
diffusion  et  on  y donne  des  rSsultats  exacts.  Nous  comparons  les  resul- 
tats  numeriques  obtenus  en  suivant  cette  theorie  avec  ceux  atteints 
par  le  biais  d'une  thgorie  approximative  en  usage.  (NC) 


ABSTRACT 


A theory  for  small  angle  multiple  scattering  of  a laser 
beam  by  aerosols  is  presented.  Unlike  in  previous  studies,  the  pre- 
sent formulation  yields  exact  results  with  which  the  scattered  radiance 
can  be  computed  to  any  desired  order  of  multiple  scattering.  Numerical 
results  based  on  the  exact  theory  are  compared  with  those  derived  from 

an  approximate  theory  currently  in  use.  (U)  — 
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1.0  INTRODUCTION 


The  performance  of  optical  communication  systems  using  laser 
sources  is  often  limited  by  the  occurrence  of  low-visibility  meteoro- 
logical conditions.  This  is  due  to  the  presence  of  a large  number  of 
aerosol  particles  in  the  atmosphere.  These  particles  scatter  as  well 
as  absorb  radiant  energy  from  an  incident  laser  beam.  As  a result, 
a well  collimated  beam  will  disperse  more  readily  in  an  atmosphere 
heavily  laden  with  aerosols  and  the  beam  intensity  will  rapidly  decrease 
correspondingly.  To  evaluate  the  operational  range  of  a given  optical 
system,  it  is  necessary  to  describe  quantitatively  the  degradation  of 
optical  signals  by  aerosols. 

It  is  well  known  that,  in  a clear  atmosphere,  the  detection  of 
laser  illumination  of  a military  target  can  be  achieved  only  when  the 
detector  is  situated  close  to  the  beam  axis.  This  is  because  the 
laser  beam  used  has  a small  angular  divergence.  To  enhance  the  off- 
axis  detectability,  it  has  been  suggested  that  aerosol  particles  may  be 
injected  into  the  atmosphere  so  that  radiant  energy  can  be  scattered 
into  a wider  angular  region. 

The  purpose  of  this  repoTt  is  to  formulate  a theory  to  calcu- 
late the  intensity  of  a laser  beam  scattered  by  aerosols.  Since  we 
are  mainly  interested  in  optically  dense  media,  our  starting  point  is 
Wentzel's  summation  method  which  includes  explicitly  multiple  scatter- 
ings of  all  orders.  Assuming  the  scattering  phase  function  to  be 
strongly  peaked  in  the  forward  direction,  Amush  and  Stotts  (Refs  1,  2) 
have  recently  studied  problems  similar  to  the  present  one.  However, 
to  arrive  at  their  analytic  results  further  approximations  were  made. 

On  the  other  hand,  these  additional  approximations  are  not  necessary 
in  our  treatment.  Thus,  for  aerosols  with  a highly  anisotropic  phase 
function,  the  present  result  is  exact. 
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In  Section  2,  Wentzel's  summation  method  is  first  discussed 
and  applied  to  the  problem  of  small  angle  multiple  scattering  of  a 
collimated  laser  beam.  It  is  also  shown  that  similar  results  can  be 
obtained  by  starting  from  the  radiative  transfer  equation.  Section  3 
contains  numerical  results  based  on  analytic  expressions  of  Section  2. 
Comparison  with  results  derived  from  applying  Arnush  - Stotts  type  of 
approximation  (Refs  1,  2)  are  also  made.  The  main  conclusions  of  this 
report  are  contained  in  Section  4. 


This  work  was  performed  at  DREV  between  March  and  November, 
1977  under  PCN  33A11,  Aerosol  Studies. 


2.0  SMALL  ANGLE  MULTIPLE  SCATTERING  OF 
A COLLIMATED  LASER  BEAM 


The  collimated  laser  beam  we  consider  can  be  regarded  as 
a monochromatic  pencil  of  light  with  a small  angular  divergence.  The 
beam  propagates  in  a medium  whose  scattering  phase  function  is  assumed 
to  be  highly  anisotropic  with  a strong  peak  in  the  forward  direction. 
This  ass'imption  is  valid,  for  example,  for  the  propagation  of  visible 
and  near  infrared  radiation  in  haze  and  fog.  Hence,  the  majority  of 
the  scattering  events  are  characterized  by  the  small  scattering  angles. 
In  our  discussion,  we  do  not  consider  the  effects  of  polarization  and, 
furthermore,  we  assume  that  the  optical  properties  of  the  medium  are 
not  affected  by  the  propagation  of  the  laser  beam. 


To  formulate  a theory  for  small  angle  multiple  scattering 
of  the  laser  beam,  two  different  approaches  will  be  followed.  On  the 
basis  of  the  summation  method  of  Wentzel  (Ref.  3),  the  intensity  of  the 
scattered  beam  can  be  written  as  an  infinite  series,  each  term  of  which 
corresponding  to  the  contribution  to  the  intensity  by  the  scattering  of 
a given  order.  This  infinite  series  can  be  summed  and  thus,  numerical 
evaluation  of  the  intensity  is  not  restricted  to  the  lower  order  terms. 
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It  is  evident  that  the  present  problem  is  a special  case  of 
the  general  question  of  the  solution  of  the  radiative  transfer  equation. 
It  is  well  known  that  even  in  the  situation  of  a plane  parallel  me- 
dium uniformly  irradiated  by  a source,  the  solution  of  the  radiative 
transfer  equation  is  non  trivial  (Ref.  4).  With  the  additional  complica- 
tion of  non  uniform  irradiation,  the  problem  becomes  more  difficult  and 
has  not  yet  been  solved  in  the  general  form.  We  will  show  that,  with 
the  small  angle  approximation,  the  radiative  transfer  equation  can  be 
solved  for  the  multiple  scattering  of  a collimated  laser  beam.  The 
solution  thus  obtained  is  identical  to  the  one  found  by  using  Wentzel's 
method . 

2.1  The  Small  Angle  Approximation 

The  small  angle  approximation  is  of  central  importance  for 
our  subsequent  discussion.  In  order  to  give  a more  precise  statement 
of  the  content  of  this  approximation,  let  us  first  introduce  a few 
quantities  and  specify  the  geometry  of  the  problem. 

We  choose  a coordinate  system  such  that  the  laser  beam  has 
its  source  at  the  origin,  and  the  beam  axis,  along  the  z-axis  (Fig.  1). 
The  intensity  I (r , s)  is  a function  of  the  position  vector  r = (x,  y,  z) 
and  a unit  vector  s specifying  a given  direction.  The  direction  s 
may  be  given  in  terms  of  the  polar  angles  (6,  <j>)  , or  the  projected  an- 
gles (A  , A ),  (Fig.  2).  It  is  clear  that  the  projected  and  the  polar 
x y 

angles  are  related  by: 

tan  <|>x  = tan  0 cos  <p  (1) 


tan 


= tan  9 sin  A 


(2) 


The  intensity  may  be  denoted  also  by  I(x,  y,  z,  9,  A)  or  I(x,  y,  z, 

Ax,  Ay) . In  the  small  angle  approximation,  it  is  assumed  that  one  can 
replace  sin  0 by  0,  and  cos  0 by  1 . In  addition,  all  the  functions,  with 
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0,  <J>  » 4V  as  arguments,  which  occur  in  I(r,  s)  are  assumed  to  decrease 

x y 

rapidly  as  6,  | d>  |,  |*  | increase  so  that,  if  g(r,  d>  , 4 ) is  such  a 

Ay  Ay 

function,  then,  the  following  approximation 
n 00 

/ (2)  g(r,  4>x , 4y)  d<J>x  d4y  = f { 2)  g(r,  4>x , 4y)  d$x  d4>y  (3) 


is  valid.  In  the  above  equation,  (2)  indicates  a 2 -fold  integration. 


Consider  a photon  with  initial  direction  of  propagation 

(4  » 4 ) at  the  origin.  It  suffers  m scatterings  without  being  ab- 

x0  y0 

sorbed  in  traversing  a distance  z,  and  the  corresponding  scattering  angles 

are  C 4Y  , 4V  ) , (4V  » 4V  ) , • • • , (<f>  , 4V  ) ; (Fig.  3).  On  emerging 
1 yl  *2  y2  m ym 


from  the  plane  z,  the  direction  of  the  photon  has  (<j>  , <f>  ) as  the  pro- 

x y 

jection  angles 


= 4 + 

T 


y y 


iu 


In  the  small  angle  approximation,  tan  4 - 4 , tan  4 - 4 , the  later- 

x y y y 

al  coordinates  (x,  y)  of  the  photon  in  the  plane  z are  given  by 


x = ^x  + (z " V 


y = + 


X tz  - IJ)  7y, 
j = l J 


where  z^  is  the  plane  in  which  the  j scattering  occurs. 
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2.2  Small  Angle  Multiple  Scattering  by  Wentzel's  Method  (Ref.  5 


Consider  a beam  of  photons  with  initial  direction  of  propa- 
gation (4>  , <j>  ) . We  choose  a plane  of  observation  z sufficiently  close 

x0  y0 

to  the  source,  so  that  only  single  scattering  needs  to  be  considered. 

Let  *(♦,♦,  x,  y,  z)  be  the  single  scattering  joint  distribution 
x y 

function  for  finding  a photon  which  has  undergone  a scattering  with 

scattering  angles  (<J>  , 6 ) and  lateral  deflections  (x,  y)  . If  the 

x y 

scattering  occurred  in  the  plane  z’,  then. 


x = (z  - z’)<f)x,  y = (z  - z')<f>y  (8] 

If  p(d  , <j>  , z)  is  the  normalized  scattering  phase  function,  the  joint 

^ y 

distribution  function  w(<j>  , <f>  , x,  y,  z)  can  be  expressed  as  (Ref.  5) 

x y 


W(-V  V,  x»  y>  z)  = 47  p(ix»  z’)  6ix  - (z  - z,Hx]  • 


6[y  - (z  - z')<t>y] 


where  is  the  scattering  coefficient. 


We  remove  now  the  restriction  that  the  plane  z is  close  to 

the  source.  The  intensity  distribution  F (4>  , <t>  , x,  y,  z)  correspond - 

x y 

ing  to  an  initial  beam  of  direction  (<j>  , 4 ) can  be  written  as  an 


infinite  series 


V TV 


F(4>x.  <fiy>  x,  y.  z)  = 'y  ^ Fj(4>x»  <ty  X,  y, 


where  Fj  is  the  contribution  from  the  j order  scattering.  The  prob- 
ability for  a photon  to  pass  through  a layer  of  thickness  z of  the 
medium  without  suffering  scattering  nor  absorption  is  given  by 


z 

‘ /° 
0 


(z')  dz' 
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where  a( z)  is  the  volume  extinction  coefficient  of  the  medium  in  the 
plane  z.  It  is  clear  that 


r z i 

'0Ox,  4y  x,  y,  z)  = IQ  exp  - f z')  dz' 

0 


«(4>  ’ ) 

0 


«WV  - 5 6(x  - ) <5  Cy  - zd>  ) (1] 

7 0 y0 

where  1^  is  a constant  proportional  to  the  photon  flux  in  the  incident 

beam  propagating  in  the  direction  Cd>  . $ ) . The  probability  of  a 

X0  y0 

photon  in  the  beam,  undergoing  one  scattering  with  scattering  angles 

(<f>  > <f>  ) and  corresponding  lateral  displacements  (x. , y-)  is 

1 yl  11 


z z 

- J cr(z')  dz'  f dz  »(TX,  Ty,  X,  y,  z') 

_ 0 Jo 


The  second  factor  in  (12)  gives  the  probability  of  one  scattering  occur- 
ring at  any  point  between  the  interval  (0,  z) , and  the  first  factor 
gives  the  probability  of  no  scattering  nor  absorption  occurring  at 

any  other  point.  From  (12)  we  can  write 

_ z z 

W V *’  Y>  Z)  = X0  exp  " f a(Z,)  dZ'  f dZ'  W(V  ♦y*  X»  y*  Z') 

0 J o 


* 1 n exP 


z z 

- f z’)  dz'  f dz'  w^x  ■ ^X  * 

0 Jo  0 


py  ‘ V x ' x°’ y ‘ y°'  z,) 


The  probability  of  having  exactly  two  scatterings  in  the  interval 
(0,z),  with  scattering  angles  and  lateral  deflections  given  by 
♦y^)*  (?  , $y  )»  (Xj , yj)  and  (x2,  y2)  respectively,  is 


A 


.V.  . 


- 
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o(z')  dz* 


L 

) dz’  . / w(?  ♦ , yv  zj  dzx 

0 

z 

/-CV  ^y2  ’ x2  ’ X2’  Z2^  dz2 


In  the  above  expression,  the  second  factor  corresponds  to  the  first 

scattering  event,  and  the  third  factor  to  the  second  scattering.  Hence, 

the  limits  of  integration  for  the  second  scattering  extends  from  z 

(where  the  first  scattering  is  supposed  to  have  taken  place)  to  z. 

As  usual,  the  expression  can  be  written  in  a more  symmetric  form  (Ref.  6) 

arJ  integrated  over  5 , i , x y to  obtain 

1 yl  1 1 


In  exP 


> <f>  , x,  y,  z)  = 

^ k y 


f °iz': 


d\d*y1dxl<i>’l  ' /dzl  "C*x,’  Vj'  xl'  >V  ZP  ■ 


dz-  w(<)>  , <|>  , x_,  y_,  z-) 

z *2  y 2 i t t 


where  <p  = $ +<p  + <f>  , <t>  = ?v  + <i>  + <J>  , x = x.  + x + xn  and 

x ^ x2  0 y yl  y2  y0 

y = yl  + y2  + V In  orc*er  to  av°id  confusion,  we  use  variables  with 
a bar  to  indicate  scattering  quantities  (scattering  angles  and  deflec- 
tions produced  by  a given  scattering) . Quantities  without  a bar,  on 
the  other  hand,  are  measured  from  the  coordinate  axes.  The  contribution 


from  j order  scattering  to  F can  be  written  similar  to  (14)  as 
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Fj (*x’  *y»  X’  y’  Z ) = *0 


00 

■ / 


(4)  d<(>x  d?y  dx2dy2 


exp  T - / , 
L 0 

j 

00 


a(z’)  dz' 


(4)  dd>  d<p  dx,dy. 
X1  yl  1 


(4)  d<J>  d*  dxi -ldyi-l 
Xj-1  yj-l  J 1 } 1 


z 

‘l  “V  Vj-  *1-  Vl-  zl>  • • • J dzj-l  “(*x,  j-  *yj  l-  *j-l' 

z _ 

y.  z.  ,)  . I dz  w(4>  , 4>  , x.,  y.,  z.)  (1! 

3-1  j-1  J x.  y’  j j 

•'r*  J -x 


where  <t>  = d> 

x x 


0 + S S 


+ E s 


X = X0  + X]  Xi  and  7 = y0  + X yi 


Let  ?(CX,  Cy,  cx,  Cy,  z)  be  the  Fourier  transform  of  F(4'x,  <Py>  x,  y,  z) 
so  that 

00 

„ f iCC<t>  + ? <t>  + £ x + c y) 

*ttx,  5y.  Cx,  V Z)  - J (4)  e x x y y x yX; 


F(<f>  # 4>  » x,  y,  z)  d<ji  dd>  dx  dy 

A / A X 


From  (10)  and  the  linearity  of  the  Fourier  transform,  we  have 


^ (^x»  ^y»  cx*  S’’  Fn ^x*  ^y*  ?x’  Cy'  ^ 
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and 


= (Cx’  5y’  Sc’  V ^ = 


From  Eqs.  (11)  and  (13) 


00 

/ 


(4)  FnWx»  x»  y»  *•)  d1>xd4>ydxdy 


(18) 


-U  (z)  i(ScSc  + W + zScSc  + zMy  5 

fc  rr  _ „ „ . “O'-2-1  X X0  X X0  X X0  y y0 

*0(Sc’  ?y’  Sc’  V Z)  " XOe  6 

(19) 


^Sc’  ^y*  Sc’  Cy*  " I~e 

z 

■f 


i(£  4>  + £ <f>  + z?  <f>  + z£  <f> 

-nQ(z)  xTx0  yyQ  ^xvxQ  Vy0 

e 


w(Cx,  5 , Cx,  Sy,  z’)  dz’ 


(20) 


where 


L 

VZ>  - / 


a(z')  dz'. 


(21) 


00 

f f i (£  ? + CJv  + c x + c y)  1 

»«x>  V V V 2,)  ■ J t4)Le  yy  y J • 


w(<t>x.  <t>y,  x,  y,  z)  d<f>xd<()ydxdy 


(22) 


Applying  the  convolution  theorem,  one  can  readily  show  that  if  ^(£), 

r\j 

g(£)  are  the  Fourier  transforms  of  the  functions  f(x)  and  g(x)  respec- 
tively, then,  the  Fourier  transform  of 


00 

/ 


f(y)  g(x  - y - yn)  dy 


a,  * ^y0 

is  f(C)  g(£)  e \ 


■MHMH i 


. 
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Consequently, 


-°n(  z)  i(EJ,  + + + t„yj 


^n^x*  V ?x*  V ^ = I0  n! 


’x*x0  yy0  x 0 yo' 


£ 

/-«.■  cy’  5x*  V z,)  dz' 


(23) 


and 


^(Cx,  ?y,  ?x,  5y,  z)  = IQ  e 


^xn(5x  + z5x}  + ^y  «y  + ZV 


eD(z)  - nQ(z) 


(24) 


with  D(z)  defined  by 
z 


fi(z) 


J w^'x’  ?y’  Cx’  Cy’  z ^ 


dz  1 


(25) 


If,  instead  of  an  unidirectional  beam,  we  have  a pencil  of  light  with 


-&(♦*  + *1  ) 

0 ^0 

a Gaussian  angular  distribution  IQ  e produced  by  the  laser 


source,  the  intensity  I(x,  y,  z,  d>  , 4>  ) is  given  by 

x y 


I(x,  y,  z,  <(>x,  4^ 
00 

/ 

-oo 

00 

/ 


00 

// 


d<t>  d(t>  F C4> x,  <f>  , X,  y,  z)  I„  e 

Xn  yn  A y u 


-6(*x  + t2v  ) 
0 y0 


(2ir) 


(4)  d5xdCyd?xdCy  e 


'0  70 

-*(Mv  + + + ^y) 


xx  y y x 


o(z)  - n0(z)  / -0(*x  + *y03 

e / (2)  d*  d*  e 


Y 


i*xft«x  + "Sc5  + ^yn^y  + z9 


X0  yo 


(26) 
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Since 


(2)  d<J>  d<j)  e 
0 y0 


i*x  (Cx  + z^x)  + i*v  (Cv  + z V -*(£  + *v  ) 

0 y0  y y 0 yo 


ir  «x  + z^2+  («y*'*'y)2 

I"  exP 46  *“ 


we  get  for  I(r,  s)  the  expression 


I(x,  y,  z,  <f>  , 6)  = 


v*  O = — ^ / (4)  dC  d£  dt  dt 

X C2tt)4  ■/  x y x y 


-i(£x*x  + Vy  + ?xx  + CyX)  n(2)  ‘ Vz) 


TT 

8 eXP 


«x  * 25x)2  * (E„  * zCy)2 


Eq.  (28)  is  the  required  result  based  on  which  the  intensity  of  a colli- 
mated laser  beam  after  undergoing  multiple  small  angle  scattering  can  be 
obtained . 

To  proceed  further,  we  make  two  simplifying  assumptions  as  in 
(Ref.  2).  Let  us  consider  the  medium  to  be  homogeneous.  It  then  follows 


« * / 


a(z?)  dz1  * az 


L 

n(z)  = j w(sx,  Cy,  cx,  c , z') 
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z _ _ _ 
f C i(5v  + + C x ♦ 

j i!'j  (,)i 


P($x>  4>y)  <5[x  - (Z  - Z ') <f>x]  <5 [y  - (Z  - z')<J>y]d$xd<i>ydxdy 


z « 

■ f iz'  f 


(2)  e 


i<(,x^x  + (z  ‘ Z'HX]  + + (z  ' z'^y^ 


■7~  p Cd>  > <P  )d 4>  d<f> 
4ir  x y x y 


(30) 


In  addition,  we  will  assume  that  the  phase  scattering  function  can  be 
approximated  by  a Gaussian 


P(*x.  *y)  ■ w e 


2 _2 
"“(♦x  + $y) 


(31) 


This  is  true  for  the  scattering  of  visible  and  near-infrared  radiation 
by  haze  or  fog.  Substituting  (31)  into  (30)  we  have 


0(z) 


4 w 

■ r-  / " dzi  f 


(2)  dd>xdd>y  e 


i*x  [ex  - c*  - z'Kx] 


2 _2 

i*y  Uy  + (Z~Z')  ?y]  ‘aWx  + V 

e 7 ’ 1 e 7 


(32) 


We  can  again  apply  the  result  of  (27)  to  evaluate  the  integrals  of  $x> 
?y  in  £2(z)  : 


5 « / •*»  ( - — 

n * 


(z  - z')cx)2  + (Cy  ♦ (z  - z')Cy)2 


4a 


| dz ' 
(33) 


Thus,  we  can  rewrite  (28)  as  follows: 
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I(x,  y,  z,  <Px,  <f>  ) 


00 

!0  Tr  _-0Z  / ,-t 

4 — e / (4)  ds< 

(2*)  P J 


d £d  £ e 


-id  t ♦ t Jp 


«*{£  ? / «p[- 

° (34) 

where  t = (?x,  Cy)  » t = (Cx  + Cy) , $ = (4>x,  <t>y)  and  rx=  (x,  y)  . In 
the  above  expression,  we  can  carry  out  analytically  either  the  inverse 
Fourier  transform  or  the  integration  over  z.  For  the  evaluation  of  the 
lowest  order  scattering  terms,  it  is  advantageous  to  carry  out  the  inverse 
Fourier  transform  analytically.  On  the  other  hand,  to  evaluate  the 
contribution  from  all  the  high-order  scattering  terms  to  I(r,  s) , the 
integration  over  z gives  a more  convenient  result. 

To  perform  the  inverse  Fourier  transform  we  expand  the  factor 


«*{£  ? / «p[-  |tMz4aI,)t|2]} 

CO  Z 

= ]C°s  /("Odv  • dzm  e*pj-  ^ 

m=n  n \ 


+ (z  - z2)c  | 


|S  + ( * - *m)C| 


z - z 
4a 


nr  I2 


f l 

S°s  if-  J (m)dzi-  • dzm  exp  " "ir  f 

m=0  0 ' L 

it.t  ZjOi.)  • Id2  Z2(m)  | 


mitl2  * 


UNCLASSIFIED 

16 


where  Z^m)  = ^ * (z  - z^), 

j = l 


a.id  Z2(m)  (z  - z.)2  (: 

j=l 

Substituting  (35)  into  (34) , we  obtain  the  following  series  expansion 

00 

I(x»  y>  2.  ♦*.  ♦y)  = ^ Im(x,  y,  z,  tx,  * ) (2 


and  Im^x’  y’  z>  ^x ’ *y-*  is  dfcfined  as 


V*.  y,  z,  * o - -S-J  TT  e'°l 

x y (2ir)4  p 


e-KE*  MrJ  y („j 

0 

|t|  Z2(m)j|.  exp  j-  ^ £ | 


oo 

m /• 

!i-  f 

m ! J 


(4)  dfdC 


(m)dZj.  . dzm  I exp 


j-M 


m | C | + 2U  Z^m)  + 


||2  + 2z(t?)  + z2|t|2 


We  interchange  the  order  of  integration  so  that 


111  ** 

y*.  y,  z,  v V =^4  -r  e'oz  TT  j • «*.  • 

00 

f (4)  itit  exp  | - r m2( 4-  . -L-)  . 2?!  ( . 4*-  ) * 


{ 
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It  can  be  shown  that 

-iitt  * ».) 


/ 


(4)  dtdc  e 


'xp  { - [l^l2  (tst  * -Jr) 


*>(^•4 0*  |t|2( 


z,(«)  2 

♦ 

4a  46 


)]| 


(2«) 

a &c£t  exp 


I 1 I*2™  , *2  ) ill2  , , z \fT  -i 

}Km)  \ ~4o  + T&  ) '*'  + 2 ^ ~4ci  + 


( — * 

1 \ 

if  I2  1 1 

\ 4a 

46/ 

K1  Jf 

(41) 


with  A(m)  given  by 


, / » . 1 W Z2(m)  . z2  \ /Zl(m)  . z Y 

( J * 4 ^ 4o  4B  ) \ 4ot  + 46  ) ' ^ 4a  46  / 


(42) 


th 


Hence,  the  contribution  to  I(r,  s)  by  the  m order  scattering  is 

I. 


m 

Im(x,  y.  z.  4>x.  V " 7T7  “T  e"°Z  5T-  C2*)5 

(2lt) 


Z 

r 


J (m)  dz, . . dz 


1"  * m A(m) 


exp 


j j_  r / z2w  + *l  ^ ni: 

j A(m)  L " \ 4a  43  / * 


2(tT  * lr  - (ir  ♦ ir  ) i^'2  ] 


(43) 


the  numerical  evaluation  of  which  involves  the  computation  of  a multiple 


th 


integration  of  m order. 


Alternatively  we  may  integrate  analytically  the  following 
integral  occurring  in  (34) 
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z 


from  which,  the  intensity  I(r,  s)  can  be  expressed  as  an  inverse  Fou- 
rier transform  given  by: 

00 


Although  equation  (43)  when  substituted  into  (38)  gives  the  same  result 
as  (45),  however,  for  numerical  computation  (43)  can  be  used  more 
conveniently  for  low-order  scattering.  We  recall  that  in  the  series 
expansion  the  number  of  integrations  to  be  carried  out  numerically  is 
equal  to  the  order  of  scattering  being  considered.  To  find  the 


T 
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contribution  from  all  the  high-order  multiple  scattering  terms,  the 
result  in  (45)  can  be  used  and  the  four-dimersional  integration  have 
to  be  evaluated  numerically.  Despite  the  fact  that  (45)  can  be 
expanded  into  an  infinite  series  with  each  term  corresponding  to  a 
specific  order  of  scattering,  the  application  of  such  an  expression  to 
calculate  the  lower  order  terms  is  not  as  convenient  as  using  (43) . 

Even  for  the  first-order  term,  a four -dimensional  integration  over  an 
infinite  domain  has  to  be  computed.  It  can  be  seen  from  the  explicit 
expression,  or  from  phyrical  considerations,  that,  for  the  scattering 
of  a tight  beam  by  an  anisotropic  medium,  the  beam  will  not  be  signifi- 
cantly broadened  by  lower  order  multiple  scattering.  Therefore,  the 

integrands  corresponding  to  these  terms  obtained  from  expanding  (45) 

— ► — ► t 

tend  to  spread  out  in  the  £,  C space  making  numerical  computation 

laborious . 


2.3  Solution  of  the  Radiative  Transfer  Equation  in  the  Small-Angle 
V Approximation 

In  the  absence  of  sources  inside  the  medium,  the  radiative 
transfer  equation  is  given  in  its  general  form  by 

S.V-+  I (r , s)  + ol (r , s)  = JJ  p(s,  s)  I(r,  s)dtt+'  (46) 

For  simplicity,  we  have  taken  the  medium  to  be  homogeneous  such  that  the 
phase  function  is  independent  of  the  position  vector  r.  If  l,  j,  k are 
unit  vectors  along  x,  y and  z axes  respectively,  the  vector  s can  be 
written  in  terms  of  the  polar  angles  (6,  <j>)  as 

s = sine  cost  i ♦ sine  sint  j + cos<f>  k 

-►  * -f-  -f-  ->» 

and  s = sine'  cost'  i + sine'  sint'  j + cose'  k 


(47) 

(48) 


Let  e be  the  angle  between  s and  s * From  the  identity 
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cos©  = cos0  cos0'  + sin0  sin0'  cos  (<f>  - $>') 


and  the  small -angle  approximation 


sin0  cos<J>  = 0cos<f>  = <p 


sind  simp  = 0sin<{>  = 


we  obtain 


cos  © = ( 1 - 


e2 

~2 1 + 90'  (cosij>  cosi)! ' + sim}>  sin<}>') 


~2 + 90'  (cosij)  cos<f>'  + siniji  sin^') 


• - -f  ( *x  - *p2  - 4-  c*y  - *;)2 


Thus,  the  phase  function  p(s,  s)  = p(e)  in  the  small-angle  approximation 
is  a function  of  the  differences  of  the  projected  angles,  i.e.  <p^  - 4>^ 
and  <p^  - <J> ^ only.  This  enables  us  to  write  the  radiative  transfer 
equation  in  the  form 


. 31  . A 31  . 31 

*x  “Sx-  ~2y~-  ~%z~  * ffI 


‘ TT-  / f^y  P(*x  - *i‘  *y  - *P  *<*■  f.  x,  ♦•)  (51) 

-CO 

Let  I(C,  c,  z)  and  p(£)  be  the  Fourier  transforms  of  I(r,  s)  and 
p(©l  respectively  such  that 


ht,  C,  z)  = fW  d«J>xd<J>ydxdy  I (r , s)  e 


i(6<f>  + rr 


-00 
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Pit)  - J (2)  d^ypC^,  *y)  e1^  (53) 

.00 

and  rx  = (x,y) . 

Substituting  (52)  and  (53)  into  (51)  we  obtain  the  following  linear 
partial  differential  equation  for  1( % , X,  z)  : 


’x  3£ 


+ ? 


y 3C, 


ht,  t.  z) « -|j.  ta.  t,  z>  * t,  z) 


4tt 


P(£)  ^(S,  c’,  z) 


(54) 


The  solution  of  the  above  equation  can  be  obtained  using  the  standard 
theory  which  we  will  describe  only  briefly.  Let  ?(£,  t,  z)  be  given  in 
the  form  of  an  implicit  solution 


Cy,  z,  1)  = Constant  (55) 

Since 

JL.  (8Hf/asx}  tay/agy) 

3^x  (a^/aT)  , 3Cy  (ay/al) 

and 

31  = -(3¥/3z) 

3z  (ai'/a!) 

these  relations,  together  with  (54)  give  the  following  equation  in  ¥: 


r ay  - 3¥  x 3¥ 

X 3«X  y + 


Pit) 


(56) 


and  is  equivalent  to  the  set  of  differential  equations 
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d£  d£ 

—5.  L 


'V 

- a 1 


On  integration  (57)  gives 


+ 5 z = C , E + c z = C 
x y x ’ ^yz  Ly 


exp/[  *r  rift  - o ] d: 


where  Cx,  Cy  and  C are  constants  of  integration, 
of  (58)  is 


V z>  h = f(Cx,  Cy,  C, 
The  fact  that  *(5X.  £ , 2,  ft  = Constant  enables 


The  general  solution 


such  that 


us  to  solve  for  C 


C ’ t,Cx-  V 


ht,  ?. 


2)  = $(f 


4, 

+ zO  exp  f i ° s 'iil t . -i  , 

q t?¥P[?+Cz-z')?J-CT| 


dz»  (62) 


Putting  z = 0 in  (62) , we  obtain 
C,  0)  = *(£) 

Therefore,  «(t  + z?)  = ?(t  + zt,  t,  0) 

Using  (63)  and  the  fact  that  the  medium  is  homogeneous, 


we  can  write 


( 

t 
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f(t  t z)  = 1(t  * zt  t,  0)  e’02  en  (64) 


where 


z 

Sw  p[  5 + t(z  - z')  j dz'  (65) 

0 


I 


It  is  evident  that  (64)  is  identical  to  (63)  obtained  from  using 
Wentzel’s  summation  method. 


3.0  NUMERICAL  RESULTS  AND  DISCUSSIONS 


In  this  section,  we  present  numerical  results  of  the  scattered 
intensity  of  a laser  beam  in  an  optically  dense  medium  using  the 
analytic  expressions  developped  in  the  previous  section.  For  these 
calculations,  we  assume  that  the  laser  beam  has  an  angular  divergence 
equal  to  approximately  0.1  mrad  such  that  the  intensity  distribution 
of  the  incident  beam  is  given  by 

2 2 

I(x,y,z  = 0,  4>x>  <)>  ) = IQe  ®^x+V^  6(x)  6 (y)  (66) 

8 2 

with  3 = 10  /rad  . The  scattering  phase  function  given  by  (31)  with 
2 2 

a = 10  /rad  is  used.  The  normalization  condition  of  the  phase 
function  we  adopt  is 


//  p (^x  » <t>y)  c4x  d(j>y  = 4 IT 


and  hence,  W = 4a. 


defined  by 


To  evaluate  the  scattered  intensity  Is(x,  y,  z,  4y) 
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Is(x,  y,  z. 


,)  = 


I (x,  y,  z,  $ , 
m J x 


♦ ) 
y 


(67) 


m=l 


where  Im(x,  y,  z,  <$>x>  is  given  by  (43),  it  is  necessary  to 
compute  numerically  the  m-fold  integrals  in  Im<  For  this  purpose, 
the  method  of  Lyness  symmetric  integration  rule  has  been  used  (Ref.  7). 
A sufficient  number  of  terms  of  I (x,  y,  z,  <p  , <j>  ) can  be  evaluated 
to  give  the  scattered  intensity  to  any  degree  of  accuracy.  Because  of 
the  rather  rapid  increase  in  the  computer  time  required  for  evaluating 
the  higher  dimensional  integrals,  our  program  was  stopped  when  a 
precision  of  5%  was  obtained.  In  Fig.  4a  to  Fig.  4f,  the  scattered 
intensity  is  plotted  against  the  lateral  displacement  for  cases  corres- 
ponding to  o=o  =10  ^cm  Fig.  5a  to  Fig.  5f  give  a similar  set  of 
S -5  -1 

plots  for  a = os  = 4 x 10  cm  . In  these  results,  the  typical  number 
of  integrations  required  for  each  point  was  between  3 and  6.  The 
smooth  curves  connecting  the  computed  points  were  computer  fits  with 
functions  possessing  two  continuous  derivatives. 


In  (Refs . 1 and  2) , an  approximative  procedure  was  introduced 
to  solve  the  radiative  transfer  equation  with  external  sources.  This 
consists  of  replacing  ft  (z)  of  (33)  by  the  following  expression: 


ftAS(z)  = a 


f i (£x  + (z-z')  gx)~  + (gy  + ( z - z ’ ) ty V 
J 4a 


dz  1 


= a 


1 /±2  -*■  2 1 -j2  3, 

z - 4^r  2 + c ? z + T-  c z -* 


(68) 


AS 

On  substituting  Q (z)  into  (34) , the  corresponding  intensity 
AS 

I (x,  y,  z,  4>  , 4>)  can  readily  be  obtained  in  a close  analytic  form: 
x y 
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AS  12a2 

I (x,  y,  z,  4>x , ^ — T exp 

y it  6a  z 
s 


a ip 

a z 
s 


12a  ,-*■  , -t,  2 

x (Y  - s z<fr) 


Thus,  the  scattered  intensity  becomes 


aq  12  a I 

ls  (x'y'z^x*V  = r~T“ 

r Trga  z 
s 


a 

a 


z 3 

s a z 
s 


" IQe"  sZ  6 (y  - z?)  (70) 

This  approximate  expression  is  supposed  to  be  valid  for  a z » 1 and 
I £ + zC  I / 2va  « 1.  However,  from  the  comparison  in  Fig.  6,  showing 
the  exact  result  with  those  obtained  using  (70)  the  inadequacy  of  the 
above  approximation  is  obvious. 


ed  intensity  versa 


s 


ttered  intensity  versu 
cm,  a = 4 X 10~5  cm 


■ 


UNCLASSIFIED 

39 

4.0  CONCLUSIONS 

We  have  derived  a theoretical  expression  from  which  one  can 
compute  exactly  the  intensity  of  a laser  beam  undergoing  small -angle 
multiple  scattering.  This  theoretical  expression  reveals  details 
whose  presence  cannot  be  obtained  from  approximate  approach.  In 
particular,  we  have  shown  the  existence  of  three,  fairly  well 
defined,  regions  around  the  axis  of  the  incident  beam.  In  the  first 
relatively  narrow  cylindrical  region,  the  intensity  of  the  scattered 
light  is  significant.  Then,  in  the  adjacent  second  region,  the 
intensity  drops  drastically.  The  second  region  is  followed  by  a third 
one  where  the  contributions  from  multiple  scattering  become  prominent. 

Application  of  the  theory  discussed  above  to  specific  types 
of  fog  is  in  progress  and  the  results  will  be  presented  in  future 
reports . 


J 
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